Efficient Magnetization Reversal with Noisy Currents 
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We propose to accelerate reversal of the ferromagnetic order parameter in spin valves by electronic 
noise. By solving the stochastic equations of motion we show that the current-induced magnetization 
switching time is drastically reduced by a modest level of externally generated current (voltage) 
noise. This also leads to a significantly lower power consumption for the switching process. 
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The dynamics of the ferromagnetic order parameter 
persists to pose a challenging problem of fundamental 
and applied nature 0. With increasing bit density of 
mass data storage devices and emergence of the magnetic 
random access memory (MRAM) concept, the speed and 
energy dissipation of the magnetization switching process 
have become important issues. In the present MRAM 
generation, magnetic bits are written by spatially ex- 
tended 0rsted magnetic fields, which sets limits to bit 
size and power consumption. An attractive alternative 
method is the current-induced magnetization switching 
predicted by theoreticians 0, El and confirmed experi- 
mentally in nano-pillar devices 0, 0, 01- In these spin 
valves, which consist basically of an electrically connected 
ferromagnet | normal | ferromagnetic metals sandwich, the 
electric current is polarized in a "fixed" layer of high 
magnetic coercivity and exerts a "spin-transfer torque" 
on the second magnetically soft ferromagnet as sketched 
in Fig. 1. Recently, time-resolved measurements of the 
current-induced magnetization dynamics have been re- 
ported ■ Advanced theoretical models for the spin and 
charge transport in magnetic devices H 0, 0] lead to 
a reasonable description of the magnetization dynamics 
within the macro-spin model, in which the magnetiza- 
tion is assumed to move rigidly under external magnetic 
field and spin-transfer torques 0, 0] . We should note 
that in larger devices evidence has been found for spin 
waves and more complicated excitations that require full 
micromagnetic simulations |l3j . 

The main obstacle that prevents wide application of 
current-induced switching is the high critical current 
needed to reverse the magnetization. There are pro- 
posals on how to reduce the critical current by sam- 
ple design |l2| and optimizing the switching process by 
a precessional switching strategy [l4j |. Since the spin- 
transfer torque vanishes for the collinear stable point of 
a spin valve, the switching time depends strongly on pro- 
cesses that induce a canting between the magnetizations 
such that the spin-transfer torque starts to kick in. This 
happens for example by increasing temperature, and is 
the basic idea of the pre-charging strategy by Devolder 
et al. 15]. More advanced strategies used in conven- 
tional magnetization switching require pulse-shaped mi- 
crowaves and rely on precise knowledge of the mag- 



netization dynamics with proper feedback. 

Electrical and magnetization noise usually degrades 
device and system performance and often efforts have to 
be undertaken to reduce it as much as possible. In nonlin- 
ear systems intentionally added noise may e.g. enhance 
the quality of signal transmission by the phenomenon of 
stochastic resonance. Noise generators find useful ap- 
plication, for example to test the response of a system 
to noise, to generate random signals for use in in en- 
cryption or to minimize the effect of quantization errors 
by a method called dithering. In this Letter, we pro- 
pose a simple method to improve the energy-efficiency of 
current-induced magnetization switching by adding noise 
to the electric circuit connected to the device. We demon- 
strate that this leads to a reversal process with increased 
switching speed and less energy dissipation. 
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FIG. 1: Schematic picture of the spin valve under consider- 
ation. The applied potential on the left/right side is ±V/2. 
The charge current I flows clockwise and G is a noise gener- 
ator. 

The magnetization dynamics in the macro-spin model 
for a small magnetic grain is described by the Landau- 
Lifshitz- Gilbert (LLG) equation [ljj 



rh = —7m x H e ff + am x m 



(1) 



where 7 is the gyromagnetic ratio and a the Gilbert 
damping constant. The magnetization direction m = 
(cos</?sin0,sin</?sin#,cos$) is parameterized by the an- 
gles 9 and 95, and we consider switching between 9 = 
(magnetization pointing up) and 9 = n (down). When a 
current bias is applied to the device in Fig. ^ the LLG 
equation for the free layer Ff rco needs to be modified to 
include the spin-transfer torque m t0 rquc = — (7/.Ms)mx 
Ifrcc x m, where Ads is the total magnetic moment of the 
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ferromagnet. The electric charge and spin currents driven 
throu gh the systems fluctuate due to thermal and shot 
noise pj| . Here we wish to investigate the effect of inten- 
tionally applied current fluctuations with spectral density 
and bandwidth controlled by an external noise generator. 
The fluctuating (spin) current through the interfaces of 
-Ffrcc creates a fluctuating spin torque that can be taken 
into account by adding a fluctuating torque to the LLG 
equation, rhfl uct = m x (h x m) , where we introduced 
an effective random field h = hg9 + h^<\> + h r r in spher- 
ical coordinates, with (hi(t)) = and (hi it) hj (i')) = 
fiijSij (t — t') for i,j G {8,(j),r}. Such a model befits 
an electromagnetic environment with a large number of 
degrees of freedom that generates Gaussian noise with 
correlation times much shorter than the magnetization 
response time. At room temperature, shot noise does 
not have to be taken into account 0. The LLG equa- 
tion in the presence of random magnetic fields leads to a 
stochastic equation for the three components of the unit 
vector m, which can be reduced to the Fokker-Planck 
equation |l7j • This allows us to access the magnetization 
reversal time as a function of applied switching current 
and current noise. 

The switchable ferromagnet is assumed to have a uni- 
axial shape anisotropy, with easy axis along z, i.e. H ff = 
Hk cos 9z. For thermal stability of the equilibrium mag- 
netization the ratio of the height of the effective energy 
barrier and the thermal energy should be high. For our 
model, this means that HkMs 3> ksT . 

We use the formalism developed in Q to find the 
charge and spin currents through the spin valve. The 
ferromagnet-normal metal contacts are assumed to have 
identical spin-dependent conductances g<* and (in 
units of the conductance quantum e 2 /h) with p = 
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We approximate the mixing 



conductance g'*, the interface-specific parameter for the 
spin-transfer torque, by (g^T _|_ <jrJ~l J /2, which is suffi- 
ciently accurate for our purposes. The spin accumulation 
in the normal metal in quasi-equilibrium can be found by 
using current conservation and spin conservation in the 
normal metal. Spin-flip in the normal metal is disre- 
garded, while that in the ferromagnet is absorbed in the 
spin-dependent conductances. Under an applied voltage 
bias V, we then find a charge current 

9 
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and a spin current leaving the free ferromagnetic layer 
for the normal metal node 



Ifrcc 
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By inserting this expression into JTT-torquc, we find for the 
spin torque contribution to the dynamics. 



^-torque 

J S = (< 
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7J5 sin 99, 

g il ) eV/(16nM s ). 



(4) 



Since this torque vanishes at 6 = 0, 7T, the absence of 
(temperature-induced) fluctuations would imply an in- 
finite switching time. Combining all contributions dis- 
cussed above, and disregarding an angle-dependence of 
the Gilbert damping due to spin pumping 19j, the LLG 
equation reduces to 



9 = —^olHk cos 9 sin 9 — 7 Jg sin 9 + hg 
ip = jHk cos 9 + h v / sin 9 . 



(5) 



In the following we focus on the effect of current fluc- 
tuations A/ (t) = I (t) — (/) with spectrum 



dte luJt (AI (t) A/(0)), 



that are generated externally on top of the thermal noise. 
The noise spectral density is taken to be a constant S w 
over a bandwidth Aw that depends on the typical time 
scales of the system: the system is not sensitive to fluc- 
tuations slower than the timescale of the magnetization 
dynamics (the inverse Larmor frequency) or faster than 
the electron transfer time defined by voltage. Let us de- 
fine the fluctuations of spin currents (where a =T,I for 
a chosen quantization axis), A/"" (t) = I a (t) — and 
the corresponding noise power 



S° 



1 



Ar (t) Ar' (*') + Ar' (t') Ai a (t) ) . (6) 



Since I (t) = I* (t) + I- 1 (t), the charge noise power can 
be written S ch = <5" T<T ■ The spin current polar- 

<vt'=T,1 

ized perpendicular to m transfers angular momentum to 
the ferromagnet almost instantaneously at the interface. 
Taking the spin quantization axis in the direction of 9, 
i.e. perpendicular to the magnetization of -Ff ree , there 
is therefore no correlation between the (transverse) spin 
currents ZJT (t) on both sides (r = L,R) of Ff rco . The 
component of the field h relevant for our discussion is hg, 
transverse to the magnetization, and can be expressed as 
he(t)= E (jh/2eM s )(lHt)-lHt))- 

r=L,R 

We describe the switching process in the Neel-Brown 
model of thermally assisted magnetization reversal (lT^ . 
We consider a statistical ensemble of particles on the 
unit sphere that represents the probability density for 
the magnetization direction. Under the influence of the 
(random) forces these particles diffuse over the sphere. 
A Fokker-Planck equation describes the evolution of the 
probability density P (9, ip, r) to find the magnetization 
vector in a certain direction [ill H?| . 
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with Ag = — jaHji cos 9 sin8 — jJss'md + ^[igg cot 8, 
Bgg = figg, A v = jHrcosQ and B vv = sin 9. 

The diffusion constant ixgg is affected by the noise power 
S w . We consider the limit of small polarization p in which 
/igg becomes independent of 9. The value of /i w is then 
equal to [igg, but is not relevant for the calculation of the 
switching time. The new diffusion parameter can then 
be written in terms of S w : 



^88 — Mthcrmal + 7T 



1 / jH 



where 



/^thermal 



2 \eM s 



Ms 



(8) 



(9) 



is the Neel-Brown diffusion constant. 

We introduce the surface probability current j (9, ip) 
defined as dW/dt = — V • j in terms of the surface prob- 
ability density W (9, <p) = P (8, ip) / sin 8. In the steady 
state 



3-0 = AgW 



1 d 



2 sin 8 88 
Since W does not depend on tp, 



BggWsinS = 0. 



(10) 



urias (a\ -iaH K sin 2 8 + 2 1 J s cos8 

W (8) oc w (8) = exp . (11) 

1^98 

In Fig. [5]the probability density is shown for some com- 
binations of current and noise power. It can be seen that 
applying a current results in an asymmetric distribution 
because the spin torque drives the magnetization to the 
second well, whereas an increased noise power broadens 
the distribution in the wells. 

For high barriers, Kramer's method can be applied to 
find the switching time of the magnet. In that limit equi- 
librium is attained separately in the regions < 8 < 8\ 
and 82 < 8 < tt (potential well 1 and 2). The effec- 
tive potentials at 8\ and 82 should be separated from the 
potential minimum by several thermal energies fc^T to 
ensure that the probability to find the magnetization in 
the middle region is low. The probability current I m in 
this middle region defined as j • 8 = I m j (2n sin 8) , is then 
small and does not disturb the equilibrium distribution 
in the two wells significantly. In the region 8\ < 8 < 82 
(see Fig. |2| then 

rii = -n 2 = 1 = -I m , (12) 

Ti t 2 

where is the probability to be in well <; — 1,2. The 
escape time from well 1, 

61 62 

ri = — / d8w(9) sin [ d£ [w (0 sin^f 1 (13) 
l^ee J J 
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FIG. 2: The probability density for the magnetization di- 
rection in the steady state as a function of the angle 9. 
The density in the absence of a current is shown for a 
noise spectral density of S w = 3.0 • 10 _20 C 2 s _1 (solid) and 
S w = 1.0-10 _19 C 2 s~ 1 (dotted). The asymmetry in the dashed 
curve is caused by a current of 0.1 mA with noise spectral 
density S w = 3.0 ■ lO" 20 ^- 1 . 



is a good estimate for the switching time. When 
2jaHx/ M-D; \°(Hk/Js\ ^ 1, we can calculate the inte- 
grals using Taylor expansions up to second order in 
and approximating sin 8 by 8 (Ref. HT 
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jaH K jaH K + 7 Js 



exp ■ 



jaH K (1 - J s /aH K y 



(14) 

We see that the spin-transfer torque reduces the effec- 
tive barrier height [ll|. while the fluctuations increase 
the effective temperature, exponentially reducing the es- 
cape time. Because of the uniaxial symmetry T2 — > T\ 
with Js — > —Js- When a clearly preferred switching 
direction is required, figg should be kept smaller than 
jJs- But when current and noise can be switched off 
immediately after a change in the magnetoresistance sig- 
nifies that switching has occured, this is no longer a con- 
straint. By combining Eq. I|14|) . © an d © we can es- 
tablish an analytical relation between the applied noise 
power, current and switching time. We use Eq. i|13|) 
to evaluate the switching time as a function of current 
and noise spectral density (see Fig. |3J). The model is 
here at room temperature with parameters typical for 
real spin valves. The ferromagnet -Ff roc is specified by an 
anisotropy field Hk = 50 mT and total magnetic moment 
Ms — 10" 17 Am 2 . The interface resistances are taken 
to be 1 fl, the polarization p of the contacts p = 0.2, 
and the damping constant a = 0.02. A typical noise 
generator power is 15 dB higher than the thermal noise 
power of a 50 resistor at 295 K. This gives a value of 
S w = 1.0 • 10~ 20 J/O. Higher noise powers, correspond- 
ing to a thermal noise at thousands of Kelvins, can be 
readily generated. 
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From Fig. [3] we see that by keeping the current fixed 
and increasing the noise, the switching time can be re- 
duced by orders of magnitude. The contours in the inset 
to the figure show how the current can be reduced by in- 
creasing the noise when the switching time is kept fixed. 
The power dissipated by the system is proportional to 

(i(tfy which increases with external noise by S^Auj/tt. 

Since the typical bandwidth Aco is of the order of 1 GHz, 
the main contribution to the power consumption comes 
from the average switching current. We observe that the 
power can be reduced by an order of magnitude. 

The calculations presented here are restricted to the 
high-barrier approximation, which sets limits to the cur- 
rents and effective temperatures. Still, the mechanism of 
noise-assisted magnetization switching should be useful 
in other regimes as well. The experimental realization of 
this mechanism might face some difficulties such as en- 
abling the high-frequency external fluctuations to reach 
the system. For instance, in the case of a large envi- 
ronmental capacitance between the leads, the impedance 
mismatch to the device might require additional mea- 
sures. Another application of this system would be the 
measurement of the noise level: by calibration of the de- 
vice, an unambiguous relation between current, switching 
time and noise power can be established. 
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FIG. 3: Switching time as a function of noise spectral density 
for several values of the current. Inset: Equal-switching-time 
contours as a function of current and applied noise spectral 
density. 

In conclusion, we propose an energy-efficient scheme 
for current-induced magnetization switching that is as- 
sisted by noise. Our approach is based on solving stochas- 
tic equations depending on the spectral density of an 
external noise source. The solution of the correspond- 
ing Fokker-Planck equation gives the dependence of the 
switching time on current and noise level. The current 
necessary to switch the magnetization can be reduced 
by applying externally generated current fluctuations. 



Without importantly complicating the device architec- 
ture the efficiency of spin-transfer torque devices can be 
improved by exponentially reduced switching times and 
an order-of-magnitude smaller power consumption. This 
could make the difference for the attractiveness of the 
current-induced switching mechanism for real-life appli- 
cations. 
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